Hint of non-standard dynamics in solar neutrino conversion 
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Motivated by the recent low-threshold measurements of the solar 8 B neutrino spectrum performed 
by Borexino, Super-Kamiokande and the Sudbury Neutrino Observatory - all now monitoring the 
transition regime between low-energy (vacuum-like) and high-energy (matter-dominated) flavor con- 
versions - we consider the role of sub-dominant dynamical terms induced by new flavor-changing 
interactions. We find that the presence of such perturbations with strength ~ W~ 1 Gf is now fa- 
vored, offering a better description of the anomalous behavior suggested by the new results, whose 
spectrum shows no sign of the typical low-energy upturn predicted by the standard MSW mecha- 
nism. Our findings, if interpreted in a 2-flavor scheme, provide a hint of such new interactions at 
the ~ 2(j level, which is rather robust with respect to 3-flavor effects possibly induced by non-zero 



PACS numbers: 14.60.Pq, 13.15. 



I. INTRODUCTION 



One of the most important achievements in our under- 
standing of neutrino properties is undoubtedly consti- 
tuted by the (indirect) proof of the existence of matter 
effects in solar neutrino flavor conversion, as predicted by 
the Mikheev-Smirnov-Wolfenstein (MSW) mechanism [J, 
2|. In fact, the explanation of the solar neutrino prob- 
lem requires a peculiar energy dependence of the electron 
neutrino survival probability (P ee ), which is elegantly 
provided by MSW transitions in adiabatic regime Q . At 
low energies (E < 1 McV) matter effects play a negligi- 
ble role and a (averaged) vacuum-like behavior emerges, 
giving rise to P ee ~ 1/2, in agreement with the de- 
terminations of the Gallium experiments SAGE and 
GALLEX/GNO 0-0. At high energies (E ~ 10 MeV) 
matter effects dominate, leading to a much stronger sig- 
nal suppression (P ee ~ 1/3), as confirmed by the ac- 
curate boron neutrino ( 8 B v) measurements performed 
by SNO |Hni and Super-Kamiokande [H-Q. This be- 
havior is indirectly corroborated by the Chlorine [l5| ex- 
periment, whose total rate gets two (inseparable) contri- 
butions from both low- and high-energy neutrinos, and 
has recently received direct confirmation by the real time 
measurements performed by Borexino [la . Ht| . able to 
monitor separately both regimes. Such a picture is now 
considered a standard framework thanks to the spectac- 
ular results achieved by the reactor experiment Kam- 
LAND [TH, which has provided the unique oppor- 
tunity to measure the solar mass-mixing parameters in 
vacuum, independently of the dynamical effects induced 
by neutrino interactions with matter. 

It is well known 0, I20l - f22| that non-standard neutrino 
interactions (NSI) , described by low-energy four-fermion 
operators O a $ ~ V a vpff of sub-weak strength e a pGF, 
especially of the flavor-changing type (a ^ /?), can pro- 
foundly modify the flavor conversion process. Interest- 
ingly, such new interactions may produce appreciable de- 
viations only in the intermediate energy region [23| de- 



scribing the transition between vacuum-like and matter- 
dominated conversions, without affecting the well es- 
tablished behaviors observed at low and high energies. 
Therefore, the accurate observation of this energy region 
is of crucial importance for pinning down potential new 
physics beyond the Standard Model. 

Intriguingly, the first low-threshold 8 B v measurements 
performed by Borexino [l7j and by the SNO low energy 
threshold analysis [ll[ (LETA), together with those pro- 
vided by the older (SK-I [lj]) and newer (SK-III Q) 
Super-Kamiokande data, point towards an anomalous be- 
havior, showing no evidence of the low-energy upturn of 
the spectrum predicted by the standard MSW mecha- 
nism. This new circumstance suggests that new interac- 
tions may be effectively at work, affecting the conversion 
of solar neutrinos in an observable way. In this com- 
munication we quantify such an expectation by showing 
that, with the inclusion of the new spectral information, 
the solar sector data (solar+KamLAND) display a non- 
negligible preference for NSI, disfavoring the standard 
MSW picture at the ~ 2a level. 



II. NOTATION 

The evolution of a two neutrino system is governed, in 
the flavor basis, by a Schroedinger-like equation 



dx 



H 



where v a is a linear combination of 
the total Hamiltonian 



and 



(1) 



and H is 
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split in the sum of the kinetic term, the standard MSW 
term [lj, [2j, and of a new NSI- induced term. The ki- 
netic term reads if kin = UKU T , where U is the real 
orthogonal 2x2 mixing matrix depending on the mixing 
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angle #12, and K is the diagonal matrix of wavcnum- 
bers ki = m3/2E (mi and E being the neutrino squared 
masses and energy respectively). In the presence of or- 
dinary matter, the standard electroweak theory predicts 
Hf y d n = diag(V, 0), where V(x) = V2G F N E (x) is the ef- 
fective potential induced by the charged current v e inter- 
action with the electrons having number density N e (x). 
For interactions with a background fcrmion / with num- 
ber density Nf(x), the new term can be expressed as [24[ 

H$£ = y/2G F N f (x)Q ^ , (3) 

where e and e' are two effective parameters which, re- 
stricting to the case of flavor-changing interactions with 
d-quarks, are related to the fundamental vectorial cou- 
plings as 

e = e£co8 023-e£ r sin023, (4) 
e' = -e^sin20 23 , (5) 

where 823 is the atmospheric mixing angle. Taking into 
account the strong upper bounds on ej^. deriving from 
the atmospheric data analysis [25|, [26| , we can safely ne- 
glect the diagonal effective coupling e'. Therefore, the 
conversion of solar neutrinos is described by the mass- 
squared splitting Am 2 = m\ — to 2 , the mixing angle 812, 
and the effective parameter e. 
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FIG. 2: Solar v e survival probability (averaged over the 8 B 
v production region) for the best fit points obtained with 
(dashed line) and without (solid line) NSI effects. 



III. NUMERICAL RESULTS 



In our analysis we have included the data from Homes- 
take [p, SAGE H and GALLEX/GNO g-Q, SK-I pj, 
the third SNO phase [l(| , and the Borexino 7 Be data lg] . 
In addition, we have incorporated the new spectral infor- 
mation provided by SNO-LETA [HI, SK-III [T|, and 
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FIG. 1: Region allowed, after marginalization of Am 2 as con- 
strained by KamLAND, separately (left panel) by solar (S) 
and KamLAND (K) data and by their combination (right 
panel). The contours refer to A\ 2 = 1 (thin solid line), 
Ax 2 = 4 (thick solid line), and A\ 2 = 9 (dotted line). 



the 8 B Borexino data (17 L We have also included the 
latest KamLAND results 19(. For the sake of precision, 
we have incorporated both standard and non-standard 
matter effects also in the KamLAND analysis. However, 
due to the low density of the Earth's crust, both have 
only a negligible impact for the range of parameters we 
are considering. Therefore, the constraints obtained from 
KamLAND do not depend on NSI. We also included NSI 
effects in the propa gat ion of solar neutrinos in the Earth 
which, as noted in [231 ], can slightly modify the regener- 
ation effect. 



In Fig. 1, we display the results of the analysis, by 
showing the allowed regions in the plane charted by 
[sin 2 0i2, e], after marginalization of Am 2 , which in prac- 
tice is fixed by KamLAND. In the first panel we show sep- 
arately the region determined by solar and KamLAND 
data. As already noticed in (27|, for non-zero (negative) 
values of e the solar data tend to prefer larger values of 
012, with an improved agreement with KamLAND. We 
also notice that the solar data taken alone tend to prefer 
non-zero NSI. As we will discuss below, this preference 
can be traced to the anomalous behavior of the solar en- 
ergy spectrum suggested by the present data. In the right 
panel we report the region allowed by the combination of 
solar and KamAND data. This plot shows that the stan- 
dard MSW case (e = 0) is disfavored almost at the 2a 
level; more exactly we find A% 2 ~ 3.6, corresponding to 
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Standard MSW + NSI 
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FIG. 3: The (ellipse-like) region determined by all data sensitive to 8 B neutrinos is reported in both panels. In the left (right) 
panel the region spanned by the theoretical model in the absence (presence) of NSI effects is superimposed. Both experimental 
and theoretical regions are plotted for A\ 2 = 1,4. 



1.9cr. The best fit is obtained 1 for e ~ —0.16, which, 
assuming maximal atmospheric mixing [30| , corresponds 
to a difference efY — e f^ = 0.23, well compatible with 
the existing experimental bounds |3lj |. 

For definiteness, we have focused on the case of in- 
teractions with d-type quarks. However, the essence of 
our results is unaltered if interactions with u-type quarks 
or electrons are considered. In all cases the absence of 
non-standard effects is disfavored at the same statistical 
level. What changes is the best-fit value of e, as can be 
expected from the proportionality of the new dynamical 
term to the fermion number density [see eq. (3)]. 



IV. DISCUSSION 

In Fig. 2 we show the solar v e survival probability (av- 
eraged over the 8 B v production region) as a function of 
the neutrino energy, for the best fit points obtained with 
or without NSI, respectively. In both cases the mass 
squared splitting is Am 2 = 7.67 x 10 _5 eV 2 , while the 
best fit values of 9 12 are slightly different, as can be un- 
derstood from the second panel of Fig. 1. It is evident 
how in the presence of NSI, the survival probability ex- 
hibits a distinctive flat behavior for E > 7 MeV, in con- 
trast with the standard MSW case which, in the same 



energy range, presents a negative slope. 

To clarify the role of the spectral information in 
favoring the flatter NSI solution, following 
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parametrize the survival probability averaged over the 
full day (daytime and nighttime 2 ) as a second order poly- 
nomial 



P ee {E) =co + Cl (E- E ) + c 2 (E - E ) 



(6) 



where Eq = 10 MeV approximately represents the energy 
where the 8 B v experiments are most sensitive. We have 
extracted the three coefficients c, from the combination 
of all the high energy experiments, taking into account 
their sensitivities. For SNO LETA this extraction is not 
feasible, as detailed information on the bin-to-bin cor- 
relations of the spectrum is not available, and we have 
directly used the coefficients' information provided by the 
collaboration itself [ll| . 

In both panels of Fig. 3 we report the (ellipse-like) re- 
gion obtained by projecting the 3-dimensional allowed 
region onto the 2-dimcnsional parameter space of the 
two coefficients (ci, C2), which encode the deviations from 
flatness. The contours correspond to the la (continuous 
curve) and 2a (dotted curve) level. The best fit coincides 
with the origin of the axes signaling a clear preference for 
a flat spectrum. In the left panel we superimpose the re- 
gion spanned in the space of these two coefficients by the 
theoretical model in the absence of NSI. From this plot 



1 Interestingly, a similar preference for flavor-changing NSI in the 
v r ) sector has been evidenced in connection with the mis- 
match of the atmospheric mass-squared splittings observed, re- 
spectively, in the neutrino and antineutrino measurements per- 
formed by minos mm. 



2 Hereafter, we neglect the small differences in the (nighttime) 
survival probability probed at the three detectors Borcxino, SNO 
and SK, induced by their different latitudes. We have checked 
that this approximation does not alter our conclusions. 
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Fig. 1). As discussed in [23, |32|, [33| this slight tension 
can be equally alleviated by the kinematical effects in- 
duced by non-zero values of the third mixing angle #13. 
It must be noted, however, that the observed flat be- 
havior of the spectrum cannot be reproduced by a non- 
zero value of 613, since this parameter induces only an 
energy-independent suppression of the 2-flavor survival 
probability [27|. Therefore, in the general three- flavor 
case we expect only a modest reduction of the statisti- 
cal preference for non-zero NSI. The robustness of the 
hint is confirmed by the 3-flavor analysis, whose results 
are shown in Fig. 4, reporting the allowed region in the 
plane [sin 2 #13, e], having marginalized away all the other 
parameters. Allowing for 6*13 > 0, the analysis still indi- 
cates a preference for the new effects at the 1.5a level. 
The future reactor searches [34| will provide a measure- 
ment of #i3 unaffected by (standard and non-standard) 
dynamical effects. Therefore, their negative (positive) 
result, will slightly favor (disfavor) the NSI hypothesis 
discussed here. 



FIG. 4: Region allowed, after marginalization of Am 2 and 
#12, by the combination of solar and KamLAND data. The 
contours refer to Ax 2 = 1 (solid line) and Ay 2 = 4 (dotted 
line). 

we learn that the standard MSW mechanism corresponds 
to a very well definite region in this plane: the variation 
of 9i2 in the range allowed by the data produces only a 
modest excursion from the best fit point, which is ob- 
tained for a negative value of c%, corresponding to the 
expected low-energy upturn of the survival probability. 
The value of the experimental A\ 2 ~ 2.1 assumed in 
this point provides a quantitative measure of the slight 
disagreement among theory and data. In the right panel 
we superimpose the region determined in the presence of 
NSI. In this case, the theoretical best fit point lies almost 
exactly at the center of the ellipse-like region, thus offer- 
ing a perfect description of the flat spectrum indicated by 
the current data. Therefore, the NSI scenario, respect to 
the standard MSW case, "gains" a negative A^ 2 ~ —2.0, 
which can be identified as the partial contribution arising 
from the better description of the energy spectrum to the 
total value A% 2 ~ 3.6 emerging from the global fit. 

The remaining A% 2 ~ — 1.6 can be traced [27| to the 
better agreement obtained in the presence of NSI, among 
the slightly different values of the mixing angle #12 deter- 
mined, respectively, by solar and KamLAND data (see 



V. CONCLUSIONS 

We have shown that the latest solar neutrino data (in 
combination with KamLAND) favor the presence of non- 
standard dynamical terms in the MSW Hamiltonian at a 
non negligible statistical level, thus hinting at new neu- 
trino interactions. We have shown how such an indica- 
tion, already present with lower statistical significance in 
the older data [U [H| , is now enhanced by the anoma- 
lous spectral behavior observed in three experiments. 
We stress that the indication we have discussed is in- 
direct, and may be confused with other possible sources 
of anomalous spectral distortions, as those induced by 
conversions into new sterile neutrino states (36[. There- 
fore, the identification of the correct sub-leading effect (if 
any) will need further corroboration not only from new 
indispensable low-energy solar neutrino measurements, 
but also from all the remaining neutrino phenomenology 
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